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ABSTRACT 

We present B VRIK images and spectroscopic observations of the z = 0.17 host galaxy 
of the compact steep-spectrum radio source MRC B1221— 423. This is a young (~ 
10 5 yr) radio source with double lobes lying well within the visible galaxy. The host 
galaxy is undergoing tidal interaction with a nearby companion, with shells, tidal tails, 
and knotty star-forming regions all visible. We analyse the images of the galaxy and 
its companion pixel-by-pixel, first using colour-magnitude diagrams, and then fitting 
stellar population models to the spectral energy distributions of each pixel. We also 
present medium-resolution spectroscopy of the system. 

The pixels separate cleanly in colour-magnitude diagrams, with pixels of different 
colours occupying distinct regions of the host galaxy and its companion. Fitting stellar 
population models to these colours, we have estimated the age of each population. We 
find three distinct groups of ages: an old population (r ~ 15 Gyr) in the outskirts of 
the host galaxy; an intermediate-age population (r ~ 300 Myr) around the nucleus 
and tidal tail, and a young population (r Ss 10 Myr) in the nucleus and blue "knots". 

The spectrum of the nucleus shows numerous strong emission lines, including [O i] 
A6300, [On] A3727, [Sn] AA6716, 6731, Ha, and [Nil] AA6548, 6583, characteristic of 
a LINER spectrum. The companion galaxy shows much narrower emission lines with 
very different line ratios, characteristic of a starburst galaxy. 

We have evidence for three distinct episodes of star formation in B1221— 423. The 
correlation of age with position suggests the two most recent episodes were triggered 
by tidal interactions with the companion galaxy. The evidence points to the AGN in 
the centre of B1221— 423 having been "caught in the act" of ignition. However, none of 
the components we have identified is as young as the radio source, implying that the 
delay between the interaction and the triggering of the AGN is at least 3 x 10 8 years. 

Key words: Galaxies: active — galaxies: interactions — galaxies: stellar content 



1 INTRODUCTION 

The onset of AGN activity in galaxies appears to be 
closely related to starbursts, possibly triggered by merg- 
ers, which lead to large increases in the amount of ma- 
terial being fed to the central black hole, thereby trig- 
gering the radio emission. Many radio galaxies show sig- 
natures of tidal interacti ons, such as tails, bridges, shells, 
and double nuclei (see iBarnes fc Hernauisd Il992l for a 
review). HST observations of powerful radio galaxies at 

* Based on observations made with ESO Telescopes at the La 
Silla Observatory, programs 70.A-0387(A) and 66.A-0288(B) 
f E-mail: H.Johnston@physics.usyd.edu.au 



z ~ 1 show disturbe d, knotty rest-frame UV emis- 
sion in the host galaxies llLong air. Best fc Rottgerindll995l : 
iBest. Loneair fc R5ttgerindll996tl 

Compact steep-spectrum (CSS) sources are bright but 
compact radio galaxies, with sizes 1-20 kpc. Unlike the ma- 
jority of powerful radio galaxies, the radio source is wholly 
contained within the envelope of the host galaxy. They have 
a steep spectral index, a < —0.5 (S v oc v a ), and a radio 
spectrum peaking below 500 MHz. CSS sources are power- 
ful (P1.4 > 10 25 W Hz" 1 ), and make up about 20% of bright 
flux-li mited low-frequ ency radio catalogues like the 3CR cat- 
alogue (| Bennett| l962t) or the Molonglo Reference Catalogue 
llLarge e tHif[^?8 l|). Th ev typically have redshifts between 
0.1-2 fsee Tou3ealll998l . for a review). The two competing 
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Figure 1. BVRIK images of MRC B1221— 423. Knots of star-formation and a tidal tail joining the galaxy to the companion are clearly 
visible. The size of each image is 22 arcsec; the separation of the galaxy and the companion is 10 arcsec (1 arcsec = 2.88 kpc assuming 
Ho = 71 kms" 1 Mpc -1 and Qm = 0.27). The orientation of the slit for the spectroscopic observations (PA=— 18?5) is shown as the line 
joining the galaxy to the companion on the B image. The labels on the V band image indicate the various regions described in the text. 



models for CSS sources are (a) that they are small in size 
because the jets are "frustrated" by interaction with dense 
surrounding gas; and (b) that they are young sources which 
will eventu ally evolve in to classical large double-lobed ra- 
dio sources (IO'Dealll998]^ Most evidence suggests that CSS 
sources are young (e.g. iFanti et alJll99fJ) . and will fade as 
they grow into large-scale radio sources. Significant luminos- 
ity evolution is required to explain the fact that the age of 
CSS sources is only 1% that of the large radio sources, yet 
20% or more of bright ce ntimetre- wavelength-selected sam- 
ples are CSS sources (e.g. | Scheueill974tlKaiser fc Alexander! 
ll997tlReadhead et alJl99rj) . The best evidence for the youth 
of CSS sources is the detection of hotspot ad vance speeds of 
0.1 to 0.3c in compact symmetric objects fsee lConwavfeOO^ 
and references therein). 

MRC B1221-423 is one of the nearest CSS sources, 
located in a host galaxy with a redshift z — 0.1706 
). The radio source has a steep spec- 



■iimnson et alJIlE 
trum (a = —0.85) and double lobes with a separation of 
only l" 5, which lie well within the envelop e of the visible 
galaxy (|Safouris. Hunstead fc Proutonl2 003j . The structure 
is suggestive of the lobe s seen in larger double sources; there 
is no evidence of a core. ISafouris et, alJ l|2003! ) estimated the 
age of the radio source to be 10 5 yr, using the tight cor- 
relation between the total kinetic power of the jet and the 



narrow-line luminosity jRawlines fc Saunderslll99ll) . At the 
distance of the galaxy, 1 arcsec corresponds to a distance of 
2.876 kpc, assuming H = 71 kms -1 Mpc -1 and fi M = 0.27. 

The host galaxy is highly disturbed, showing signs of 
tidal interaction with at least one close companion. An 
archival image taken using the Anglo- Australian Telescope 
in l"4 seeing shows faint concentric arcs to the north-east 
JSafouris et al.ll2003l) . evidence of past tidal interaction. 

We have obtained multi-colour images of 
MRC B1221— 423 in order to study the change in properties 
of the stellar population across the galaxy. The use of 
spatially-resolved colours to study differences in stellar 
populations across a galaxy has been shown to be an 
effective way to study age distributions, localised extinction 
and star formation histor y in galaxies. The technique 
was pioneered by IBothunI 1^986), who used a spatially 
resolved (B — R) vs. B diagram to examine the star forma- 
tion h istory of NGC4449. More recently, lAbraham et alJ 
( 1999) used pixel colour-magnitude diagrams to investigate 
the evolutio nary history of ga l axies in the Hubble Deep 
Field, while iKong et all <200OT . Ide Griis et all d200St and 
iKassin et al.l J2003h used such diagrams to study various 
nearby galaxies. 

In this paper, we present and analyse five-colour images 
and optical spectroscopic observations of MRC B1221— 423. 
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In we describe the observations and their reduction. 
In q3.1l we analyse the images, first by means of colour- 
magnitude diagrams, and then by fitting spectral energy 
distributions to each pixel. In § 13. 21 we examine the surface- 
brightness profile of the galaxy. In >!3.3l we discuss the spec- 
troscopic results, and then discuss the implications for the 
history of the galaxy in iJH 



2 OBSERVATIONS 

2.1 Imaging 

Four-colour optical images were taken on 2002 February 25 
using SUSI2 on the New Technology Telescope (NTT) of 
the Europea n Southern O bservatory with Bessell B, V, R 
and /filters (Bcsscll 1990). The spatial scale on the detector 
(two 2k x4k EEV CCDs) is 0.161 arcsec/pixel. The condi- 
tions were photometric, and details of the observations are 
give n in Ta bic U The equatorial standard star field PG0942 
from lLandoltl (|W92) was observed to derive the photometric 
transformations onto the Johnson-Cousins system. Standard 
reductions were performed using the iraf software suite. 

An infrared image was obtained on 2002 December 23 
through the Ks filter using the Sofl instrument on the NTT. 
The field of view using the Large Field objective was 4.9 
arcmin, with a spatial scale of 0.292 arcsec/pixel. The auto- 
jitter mode was used to obtain fifteen separate images, with 
each image being an average of nine 7-s images; thus the 
total exposure time on source was 945 s. The telescope was 
moved by a random amount within a 100 arcsec box be- 
tween images. Reduction of the images was done using the 
xdimsum 1 package in iraf. The data were flattened using 
a sky flat created by masking out sources and combining 
the masked images; the final mosaic was then created usingj 
xmosaic. The Persson standard star SJ 9154 iPersson et alJ 
1998) was used to derive the photometric transformation. 

Finally, all five images were registered onto the same 
pixel scale, by fitting the positions of stars common to each 
image with linear transformations, including terms for trans- 
lation and scaling. The K s-band image was re-binned to the 
same scale, and a term for rotation also included. The im- 
ages were convolved with gaussians to match the image size 
on the worst-seeing image (the Jfs-band image, with mea- 
sured FWHM of 1.0 arcsec), and then binned 3 x 3, so the 
pixel scale of the final images is 0.483 arcsec/pixel, which 
critically samples the seeing disk. The registration of the 
images is accurate to about 0.2 pixels, as measured from 
the positions of field stars in the final images. 

The final set of images is shown in Figure^ The B-band 
samples rest-frame wavelengths in the range 3340-4180A, or 
essentially rest-frame [/-band. 

2.2 Spectroscopy 

MRC B 122 1—423 was observed on four occasions, in 2001 
April and 2003 August, using the Double Beam Spectro- 

1 DIMSUM is the Deep Infrared Mosaicing Software pack- 
age developed by Peter Eisenhardt, Mark Dickinson, Adam 
Stanford, and John Ward, and is available via ftp from 
ftp : //iraf . noao . edu/iraf /extern/ xdimsum/ 



Table 1. Journal of observations of MRC B1221-423. The imag- 
ing observations were taken with the New Technology Telescope 
at the European Southern Observatory; the spectroscopic obser- 
vations were taken using the Double Beam Spectrograph on the 
ANU 2.3m telescope at Siding Spring Observatory. The columns 
show the date of the observation, the wavelength range and the 
resolution in each of the spectrograph arms, and the exposure 
time. The seeing for the spectroscopic observations was 1''2— 1'/5. 



Imaging observations 



UT Date 


Instrument 


Filter 


^cxp 

00 


Seeing 
(") 


2002 Feb 25 


SUSI2 


B 


600 


0.9 






V 


600 


0.7 






R 


400 


0.6 






1 


900 


0.8 


2002 Dec 23 


Sofl 


Ks 


945 


1.0 





Spectroscopic observations 






UT Date 


A ranj 


;e (A) 


FWHM (A) 


tcxp 




blue 


red 


blue 


red 


(s) 


2001 Apr 16 


3180-5800 


5880-7510 


4.6 


2.3 


7200 


2001 Apr 17 


3180-5800 


6610-8550 


4.6 


2.3 


5400 


2001 Apr 18 


3180-5800 


6610-8550 


4.6 


2.3 


5400 


2003 Aug 3 


3180-10370 




9.7 




2400 



graph on the ANU 2.3-m telescope at Siding Spring Obser- 
vatory. The detectors were two SITe 1752 x 532 CCDs with 
15-pm pixels. A complete log of observations is given in Ta- 
bic For the observations in 2001, a dichroic filter with a 
cross-over wavelength of 5500 A was used to split the light 
into the two arms of the spectrograph. The observation in 
2003 August was made with a low resolution grating and a 
plane mirror in place of the dichroic, so that all the light 
was sent to the blue arm of the spectrograph in order to 
obtain complete spectral coverage. Pairs of 1200 s or 1800 s 
exposures were bracketed with CuAr arc-lamp exposures. A 
slit of width 1.5 or 2 arcsec was used, oriented at a position 
angle of —18? 5 so as to include light from the companion on 
the slit. On 2001 April 16, the position angle was —16°. The 
spatial scale on the detector was O'.'91/pixel, corresponding 
to a linear scale of 2.6 kpc per pixel. 

The iraf software suite was used to remove the bias 
and pixel-to-pixel gain variations from each frame. As we 
had multiple consecutive observations of the same object, 
cosmic ray events were removed using the technique de- 
scribed bv ICrokel 1^995), as implemented in FIGARO. The 
spectra were straightened in FIGARO so that the dispersion 
ran exactly along rows of the image, then a two-dimensional 
wavelength fit was performed to the arc images by fitting a 
third-order polynomial to the arc wavelengths as a function 
of pixel number, for each row of the image. These wavelength 
solutions were copied to the object images, interpolating be- 
tween the bracketing arc exposures, and the data were re- 
binned so the wavelength-pixel relation was linear and uni- 
form across the image. The sky background was subtracted 
from each image, by fitting to the sky on either side of the 
galaxy, chosen to be well outside the wings of the galaxy pro- 
file. The spectra were corrected for atmospheric extinction, 
the telluric absorption features were removed by comparing 
with the spectrum of a smooth-spectrum standard taken at 
similar airmass, and the spectra were flux-calibrated. The 
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Figure 2. (a) V vs. V — I colour-magnitude diagram for individual pixels in B1221— 423. The units are magnitudes per pixel, where 
each pixel is 0.483 arcsec on a side. All pixels brighter than V = 25 are shown as dots. Three distinct "plumes" of pixels are seen; these 
correspond to the three regions in (b). The sky brightness was V = 22.8 mag pixel -1 . The dashed line shows the relation (V — I) + V = 
24.6 mag, which we took as the boundary delimiting the regions, (b) Image showing the three regions identified in the colour-magnitude 
diagram, overlain with contours of the B-band image (Fig. 0. Light grey pixels are from region 1, dark grey from region 2, and black 
from region 3. The letters refer to model spectra shown in Fig.|I] representing pixels in the nucleus (N), companion (C), knot (K) and a 
red region to the north-east of the nucleus (R). 



nights were not all photometric, so the flux calibration can 
only be regarded as approximate. 



3 ANALYSIS 

3.1 Modelling pixel colours 

3.1.1 Colour-magnitude diagrams 

Having created a uniform set of images, with the same scale, 
orientation and seeing, we can now use this dataset to study 
the variation in the stellar population across the face of the 
galaxy. 

The pixel colours were first corrected for 
Galactic extinction, using the extinction map of 
ISchleeel. Finkbeiner fc Davis! il998T) . which indicates 
A v = 0.331 mag at the position of MRC B1221-423. 
The other wavebands were c o rrecte d using the redden- 
ing law of iRieke fc Lebofskvl <1985ft . The (V - I) vs. 
V colour-magnitude diagram for individual pixels of 
MRC B1221-423 is shown in Figure ^fa). Three distinct 
"plumes" of pixels are apparent in this diagram. We have 
mapped these three features, bounded on the lower edge by 
pixels satisfying the relation (V — J) + V < 24.6 mag, onto 
the image. This defines three distinct regions, shown in 
Figure The features in the colour- magnitude diagram 
are clearly separate in the image. Pixels in the bluest 
region, Region 1, correspond to the interacting companion, 
with a hint of a trail of pixels back up the tidal tail (light 
grey pixels in Figure Hjfft)). The reddest pixels, Region 3, 
come from the centre of the galaxy (black pixels); while 
the intermediate plume of pixels, Region 2, corresponds 



to regions at larger radial distance from the centre of the 
galaxy, including the "knots" in the B-band image (dark 
grey pixels). 

The lower boundary line in Figure l2Ta) is arbitrary, 
and was chosen to show clearly the spatial separation of 
pixels of similar colours. Other features, such as the tidal 
tail, also occupy distinct regions of the colour-magnitude 
diagram, though the separation becomes blurred as the pixel 
flux density drops significantly below the sky level (22.8 mag 
per pixel in the V band). The same separation is seen for 
other colour-magnitude combinations, but is strongest for 
widely-separated bands. 

Thus the colour-magnitude diagram separates pixels 
into distinct groups, with different regions of the diagram 
occupying different areas of the galaxy and its interacting 
companion. The most obvious interpretation of this effect is 
that the different groups represent different stellar popula- 
tions; in this case, the three groups we define here probably 
represent the products of three distinct episodes of star for- 
mation. 

3.1.2 Fitting spectral energy distributions 

A colour-magnitude diagram does not use all the informa- 
tion we have available. With five colours, we can go further 
in analysing the stellar population: by fitting spectral mod- 
els to the colours of each pixel in turn, we can determine the 
age and reddening of the stellar population across the face 
of the galaxy. 

We used a purpose-written program to model the spec- 
tral energy distrib ution (SEP) of pixels, loosely bas ed on 
the hyperz code of lBolzonella. Miralles fc Pelld (l200Cil . The 
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Figure 3. Histogram of ages for each pixel, derived from the 
best-fit model. The top panel shows ages derived from a model 
with exponentially decaying star formation with a decay time of 
t — 15 Gyr. The pixels are grouped into three distinct peaks, 
with ages of roughly 13 Gyr, 500 Myr and 10 Myr. The bottom 
panel shows a fit (excluding the pixels in the 15 Gyr peak) using 
a model with r = 100 Myr. The two young populations are still 
present, with ages around 210 Myr and 10 Myr. 

method is as follows. An input stellar population model is 
chosen, and the redshift is fixed at the known redshift of the 
galaxy, z = 0.1706. A set of model spectral energy distribu- 
tions is constructed by folding the input spectrum through 
the BVRIK response functions, resulting in a grid of model 
fluxes, with varying age and reddening along the axes, to 
be compared to the observed B VRIK fluxes for each pixel. 
The difference between the sum of the model fluxes and the 
observed fluxes in each bandpass is computed, and the pa- 
rameters of the best fit (age a nd reddening) record ed. 

Motivated by the work of lKassin et alJ <l2003t) in mod- 
elling the stellar populations in the merging galaxies known 
as the "Antennae", we used as our initial model an old 
star-forming disk model with exponentially decaying star 
formation on a timescale of r = 15 Gyr. The mo del was 
from the GISSEL96 library feruzual fc Charlotlll993l) . with a 
Salpeter initial mass function with mass limits of 0.1 and 
100 M q and solar metallicity, using the iGunn fc StrvkeJ 
(1983) stellar spectral atlas. A histogram of the fitted ages 
shows three distinct peaks, with ages of 13 Gyr, 500 Myr, 
and 10 Myr (Figure^. We can refine this by fitting a model 
with r = 100 Myr to the "young" pixels, excluding those 
with ages above 5 Gyr (Figure lower panel). The young 
pixels still show two distinct age groupings, with ages of 
210 Myr and 10 Myr and younger. 

Examples of the fitted spectral energy distributions are 
shown in Figure [I] The observed spectral energy distribu- 
tions for the four points indicated with crosses in Figurel2Kfe) 
are plotted, together with the best-fit model spectra from 
the dual-age fit described above. Point R is in the "old" re- 
gion, as indicated by the fit with the r = 15 Gyr model; the 
other three points were fitted using the r = 100 Myr model. 
The fitted parameters are shown in Table [5] 

We can examine the location of the populations of dif- 
ferent ages by plotting an age map of the galaxy. This is 
shown in Figure [K] The ages of pixels clearly correlate with 
morphological features of the galaxy (shown as contours for 
comparison). The old population is found through most of 
the galaxy, and clearly represents the underlying stellar pop- 
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Figure 4. Spectral models fitted to observed spectral energy dis- 
tributions, for the four points denoted with a letter in Figurel2Yfe ). 
and an exponentially-decaying star-formation model fitted, as de- 
scribed in the text. The squares represent the pixel fluxes to which 
the model SEDs are fitted. The horizontal bars represent the 
widths of the filters, the vertical bars the photometric errors on 
the pixel fluxes. The spectra are offset in the vertical direction for 
clarity. The pixels shown are labelled with letters in Figure l2Yfe): 
they represent pixels in the companion (C), one of the blue knots 
to the west of the nucleus (K), the nucleus (N), and a red region 
to the north-east of the nucleus (R). The parameters of the fit 
are shown in Table |21 



Table 2. Parameters of the fits to the spectral energy distri- 
butions for the four selected points shown in Figure l2Yb ): the 
fits are shown in Figure [1] We fitted an exponentially-decaying 
star formation model with a timescale of either r = 15 Gyr or 
r = 100 Myr, as shown in the second column. The best-fit age 
and reddening for each set of pixel fluxes is shown in columns 3 
and 4. 



Region 


Model 


Age 


A v 






(Myr) 


(mag) 


Companion 


100 Myr 


290 


0.1 


Knot 


100 Myr 


28 


1.5 


Nucleus 


100 Myr 


128 


2.0 


Red region 


15 Gyr 


15500 


0.6 



ulation of MRC B1221— 423. The intermediate age popula- 
tion, with ages around 500 Myr, is found in the companion 
and along the tidal tail joining the galaxy to the companion, 
while the youngest population, with ages ~ 10 Myr, is found 
in the centre of MRC B1221— 423 and in the knots surround- 
ing it, with a small group of young points in the companion 
galaxy. In the B-band, the old population is contributing 
30% of the light, the intermediate-age population 46%, and 
the young population 24%. 

The fit is quite poor in the nuclear region (point N), 
especially for the V and R bands. Examining the grid of 
models, there are two solutions which are almost equally 
acceptable for points in this region. The model shown in 
Figure 2] has an age of 128 Myr and a reddening Ay — 
2.0 mag (xS = 1-7); a second solution, with a reduced x 2 
only marginally higher (xi = 1-9): nas an a g e °f 9 Myr 
and a reddening Av = 2.85 mag. Since BVRI filters overlap 
and hence are not independent of each other, the Xv values 
cannot be taken as a quantitative measure of the goodness of 
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Figure 5. Age image of MRC B1221-423, created by fitting mod- 
els with exponentially decaying star formation, with r = 15 Gyr 
and t = 100 Myr (see text for details). Pixels in the three peaks 
of the age histogram (Fig. |3J are shown as different colours. The 
oldest pixels (ages > 5 Gyr) are shown in black, pixels with ages 
between 80 Myr and 5 Gyr are shown as grey, and the young pix- 
els (ages < 80 Myr) are shown as white. Contours of the B-band 
image are overlain for comparison. The ages of pixels clearly cor- 
relate with morphological features of the galaxy. 



fit; however, the relative values indicate the second solution 
must be almost as good as the first. 

The split between these two solutions can be seen in 
Figure |5J where some pixels in the nuclear region are white 
(young, with ages < 80 Myr), while others are grey (inter- 
mediate age, with ages between 80 Myr and 5 Gyr), despite 
having similar colours (Figure |2J • This may well reflect a 
real mix of ages in the nuclear region, so that a single age 
population does not well describe the observed colours. Al- 
ternately, it could represent a spread in metallicity within 
the nuclear region, since changes in both age and metallicity 
can lead to similar colours. 

The spectrum of the galaxy shows many emission lines 
(§ 13.311 : however, as they are comparatively weak (contribut- 
ing at most 7% of the flux) , they do not seriously distort our 
analysis. 

We also derive the extinction across the galaxy from 
the same model fits (Fig. [(J. Most of the region has Ay < 1, 
but the centre of the galaxy and the blue regions to the west 
show evidence for large amounts of dust, Av ~ 1-2. This 
means that the region with the youngest population - the 
nucleus of the galaxy - requires the most internal reddening 
to fit the observed colours. 



3.2 Surface photometry 

Since B1221— 423 is a powerful radio source, we would ex- 
pect its host galaxy to be an elliptical galaxy; yet the optical 
image shows a galaxy which looks like a spiral. Analysis of 





A v (mag) 

Figure 6. Extinction image of MRC B1221-423, from the same 
model fits as were used in Fig. |3 Black pixels correspond to heav- 
ily reddened regions, Ay > 3. Contours of the B-band image are 
overlain for comparison. 



the light profile of the galaxy seemed to be the best way to 
identify the shape of the underlying stellar population. 

We fit elliptical isophotes to the (unbinned) galaxy im- 
ages, using the task ellipse in iraf. We then attempted to 
fit either de Vaucouleurs r 1 ^ 4 -law profiles, appropriate for an 
elliptical bulge, or a combined bulge plus exponential disk 
profile, as expected for a spiral galaxy. 

The results of the analysis are shown in Figure |7| The 
-fsTs-band image is well described by an r 1//4 -law profile, but 
the V and i?-band profiles are not. This implies that the 
host galaxy does have a bulge, but it is being masked in 
bluer colours by the younger populations. We constrained 
the V and i?-band profiles to have a bulge of the same size 
as the Ks-haxiA, with an effective radius R e — 92 pixels 
(= 14.8 arcsec). A de Vaucouleurs profile with this radius 
plus an exponential disk is a good fit to both the R and V- 
band light distributions (the _R-band light profile is shown in 
Fig. . This bulge radius corresponds to a half-light radius 
of 42 kpc. Correcting for the extinction of Av = 1.9 mag as 
derived in the SED fitting, and inclu d ing K - correct ions for 
the redshift z = 0.17 from iPoggiantil (|1997D we can derive 
the intrinsic luminosity and colour of the bulge. The bulge 
has V—K — 3.1, and a luminosity in the V-band of 10 12 Lq. 



3.3 Spectroscopy 

We constructed a single spectrum (Fig. |S| by summing the 
spectral images from 2001, after aligning the spectra; we 
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Figure 7. iCs (filled squares) and ij-band (open circles) luminos- 
ity profiles for MRC B1221-423, plotted against r 1 / 4 . The solid 
line shows the best-fitting de Vaucouleurs profile to the ifs-band 
data between 4 and 20 pixels (1.4 < r 1 / 4 < 2.1). The dashed line 
shows a bulge+disk model fit to the i?-band data, constraining the 
bulge to have the same radius as in the if s-band, R e = 92 pixels. 



could only use the data from 2001 April 17 and 18 for the 
red spectrum, since we had a different wavelength setup for 
the first night (Table [TJ. 

Emission is detected along the slit from the galaxy and 
companion. The heliocentric redshift of the emission lines, 
determined from a combined fit to the red side of the nu- 
clear spectrum, is z — 0.17068 ± 0.00002. The lines show 
significant rotation, with an amplitude of 90 kms~ along 
this axis. The companion galaxy is found at a velocity of 
— 105 ± 20 kms -1 with respect to the nucleus. The line ra- 
tios can also clearly be seen to change along the slit, re- 
versing from the nucleus to the companion, with the nuclear 
spectrum showing Ha weaker than [Nil] (log [Nil] /Ha = 
+0.25), while in the companion [Nil] is weaker than Ha 
(log [Nn]/Ha = -0.35). 

We extracted the spectrum of the nuclear region by se- 
lecting the central two rows, corresponding to a linear size 
of 5.2 kpc. The resulting spectra are shown in Figure|S] The 
line properties were measured usi ng the specfit package 
implemented in iraf llKrisslll994h . Gaussian profiles were 
fitted to the lines, with a single velocity shift for all lines. 
The line fluxes and equivalent widths are given in Table [3] 

A spectrum was extracted for the companion galaxy, by 
summing together seven rows at the position of the galaxy. 
This spectrum is also shown in Figure |5] This spectrum is 
clearly very different from the spectrum of the nuclear re- 
gion, with much narrower lines and very different line ratios. 

The nuclear spectrum shows the strong [O i] A6300, 
[On] A3727 and [Sn] AA6716, 6731 emission lines char- 
acteristic of the class of active galactic nuclei known as 
"low-ionisation nuclear emission- line regions", or LINERs 
jHeckmanlll980l) . We can further investigate the classifica- 
tion of the different r egions of the galaxy by using the clas- 
sification diagrams of IVeilleux fc Osterbrockl lll987l) . 

From the ratio of the intensities of the Ha and H/3 lines, 
we can estimate the extinction towards the line-emitting re- 
gions. For the nucleus of the galaxy, the observed ratio is 10, 
compared with the expected ratio of 2.86 for a temperature 
T = 10 4 K, which implies an extinction Av = 3.7 mag. For 



the companion, the ratio Ha/H/3 = 7.6, which implies an 
extinction Av = 2.9 mag. These are both higher than the 
extinctions we derived from modelling the spectral energy 
distributions (Table where we found Ay = 1.9 mag for 
the centre of the galaxy, and Av = 0.5 mag for the com- 
panion (plus 0.33 mag of local extinction). However, since 
the companion has both Balmer emission and absorption, 
the ratio of Ha to H/3 may not be a reliable measure of the 
extinction. 

The nucleus of MRC B1221-423 has the following line 
ratios, corrected for the reddening as determined from the 
Balmer decrement: [O l]/Ha = 0.3, [N ll]/Ha = 1.9, [S ll]/Ha 
= 1.0, an d [O hi] /H/3 = 6.5; these rati os match a LINER 
spectrum JVeilleux fc Osterbrockl Il987l) . The ratio of the 
[O ill] to H/3 flux is higher than typical for a LINER spec- 
trum, but since the [O in] flux is not well determined (being 
seen only in the low-resolution spectrum from 2003 August) 
we give more weight to the high [Nll]/Ha ratio. 

For the companion galaxy, the corrected line ratios 
are: [Ol]/Ha = 0.07, [Nll]/Ha = 0.5, [Sll]/Ha = 0.4, and 
[O iii]/H/3 < —0.05 (based on a marginal detection of H/3 
and an upper limit to [Oiii]), which match the spectrum of 
a starburst galaxy. 

The limited signal-to-noise ratio of our spectra preclude 
a more detailed comparison with model population spec- 
tra, and do not allow us to make direct comparisons with 
the models derived from the broad-band images. However, 
there are a few indicators of the underlying stellar popula- 
tion which we can look at. 

The A4000-A break is a useful age-indicator for inte- 
grated stellar populations. Older populations produce a dis- 
continuity at 4000A due to the increase in stellar opacity 
produced by metal lines shortward of 4000 A in late-type 
stars. We measured the A4000-A break index by compar- 
ing the mean flux density in two windows above and below 
4000 A, taking the ratio of the mean flux between 4050 
and 4250 A w ith the mean flux between 3750 and 3950 A 
jBruzuall983t) . For the nuclear spectrum, this index is 1.37. 
This small break index confirms the exist ence of recent star 
forma tion in the nucleus of the galaxy (e.g. Kauf fmann et all 
120031) . 

Since the galaxy contains a powerful radio source, some 
fraction of the excess blue light may arise from the AGN. We 
detect the G band in the nuclear spectrum, with an equiv- 
alent width of 4.lA (Table [3J, which is consis tent with the 
value s seen in normal elliptical galaxies (e.g. iTraeer et alJ 
1998); this means that we are definitely detecting stars. The 
detection of H<5 absorption, with an equivalent width of 1.2A 
means that at least some of the blue excess is being con- 
tributed by intermediate-age stars. We conclude that it is 
unlikely that the AGN contributes more than about 50% of 
the light. 



4 CONCLUSIONS 

We have shown that the host galaxy of MRC B1221-423 
and its interacting companion can be resolved into different 
stellar populations. We have achieved this in two different 
ways, using colour-magnitude diagrams, and by using mod- 
els for the spectral energy distribution of individual pixels. 
Both these methods separate the pixels into distinct regions 
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Figure 8. Combined spectrum of the nucleus (top spectrum) and the companion (bottom spectrum), with emission and absorption lines 
identified. The nuclear spectrum has been offset in the y-direction by 2 X 10 -16 ergs -1 cm -2 A 1 for clarity. The region between the 
dashed lines is the lower resolution spectrum from 2003 August, which has much lower signal to noise; the data on either side are the 
high-quality summed spectra taken in 2001 April from the blue and red arms of the spectrograph. 

Table 3. Properties of detected emission lines and absorption lines in the spectra of the nucleus and companion galaxy (Fig. [8). The 
columns show the line with its rest wavelength, its equivalent width, integrated flux, and FWHM for the nuclear and companion spectra 
respectively. Line properties were measured by fitting Gaussian profiles to the spectra; the values in the table are the properties of the 
fitted Gaussians. The quantities with no errors (shown in italics) are linked to the corresponding quantities of another line, e.g. the flux 
ratio of the [O II] AA3727, 3729 lines was fixed at 4.3, and their FWHMs constrained to be the same. Negative values denote emission 
lines. The [Om] AA 5007, 4959 lines (marked with *) were present only in the low-resolution spectrum of 2003 August. Lines which were 
not detected are indicated by a 3-<r upper limit in the equivalent width column, calculated by EWij m = 3\/2 AA/SNR, where AA is the 
dispersion (A/pix) and SNR is the signal-to-noise ratio of the spectrum at the location of the line; limits on doublets are shown as a 
single entry. 
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[On] 3727 


-40.9 ± 1.1 


-7.3 ± 0.2 


685 ± 20 


-7.5 ± 1.4 


-2.2 ± 0.4 


500 ± 130 


[On] 3729 


-9.4 


-1.7 


685 


-1.8 


-0.5 


500 


H/3 4861 


-3.1 ± 0.9 


-1.3 ± 0.4 


360 ± 120 


-1.9 ± 1.0 


-1.0 ± 0.6 


190 ± 150 


[Om] 4959* 


-4.3 


-3 


800 


< 6 






[Om] 5007* 


-12.9 ± 1.6 


-9 ± 1 


800 ± 110 






[Oi] 6300 


-4.9 ± 0.5 


-3.4 ± 0.4 


560 ± 60 


< 0.7 






[Nil] 6548 


-11.1 


-8.5 


580 


-3.3 


-1.2 


170 


Ha 6563 


-21.0 ± 0.9 


-13.8 ± 0.4 


480 ± 20 


-18.0 ± 0.7 


-7.6 ± 0.3 


145 ± 6 


[Nn] 6583 


-33.2 ± 0.6 
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580 ± 15 
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170 ± 15 


[S n] 6716 


-11.2 ± 0.5 


-8.3 ± 0.4 


560 ± 30 


-5.2 ± 0.9 


-1.9 ± 0.3 


200 ± 30 


[Sn] 6730 


-8.0 ± 0.5 


-5.9 ± 0.4 


560 


-3.8 ± 0.7 


-1.4 ± 0.3 


200 


Can H 3933 


5.8 ± 0.7 


1.3 ± 0.1 


1080 ± 130 


3.5 ± 1.9 


1.4 ± 0.8 


790 ± 420 


Can K 3968 


5.8 


1.4 


1080 
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1.9 


790 


HC 3889 


< 1.1 






3.2 ± 1.8 


1.3 ± 0.7 


860 


He 3970 


< 0.9 






< 1.5 






H<5 4101 


1.2 ± 0.8 


0.3 ± 0.2 


1080 


5.5 ± 1.8 


2.5 ± 0.8 


860 ± 180 


G band 4304 


4.1 ± 0.7 


1.2 ± 0.2 


1080 


< 0.7 






Na Dl 5889 


1.5 ± 0.3 


1.0 ± 0.2 


620 ± 180 


< 0.8 






Na D2 5895 


0.8 


0.5 


620 
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with different coiours and/or ages. We interpret this as in- 
dicating that the gaiaxy and its companion have undergone 
severai episodes of star formation, possibly associated with 
the tidal interaction between the two. 

The ages we deduce for the stellar population fall into 
three distinct groups. An old population, with age t ~ 
15 Gyr, occupies the outskirts of the host galaxy. An inter- 
mediate age population, with t ~ 300 Myr, is found near the 
nucleus and around the tidal tail joining the two galaxies, 
while the youngest population, with ages less than 10 Myr, 
is concentrated in the nucleus and the blue regions to the 
west, as well as in the companion galaxy. 

The location of these populations is not identical to the 
separation suggested by the regions in the colour-magnitude 
diagram (Figure |5J. This may be due to the effects of dust, 
which our SED modelling suggests is concentrated in the 
nucleus of the host galaxy. 

The relationship of the different stellar populations 
with the underlying powerful radio source is much less 
clear. The age of the source, estimated from the power of 
the radio jets and to tal stored energy in the syn chrotron 
plasma of the lobes ([R ailings & Saunders! Il99lh . is only 
10 s yr ijSafouris et al]l2003h . This is much younger than 
the youngest population we detected, with an age of about 
10 Myr (Section I3.1.2fl . and suggests there must be a sub- 
stantial time delay between the most recent burst of star 
formation and the triggering of the radio source. 

A possible scenario to explain the three component pop- 
ulations is as follows: the tidal interaction began 300 Myr 
ago, which triggered star formation in both the host galaxy 
and the companion. We see the results of this episode of star 
formation as the intermediate-age population. Gas is driven 
down into the nucleus on timescales of ~ 10 s y, where it 
both cools to form stars, which we see as the 10 Myr-old 
population, and triggers activity in the central supermassive 
black hole. Delays of this order ar e predicted in theoretical 
calculations, (e.g. iLin et al]ll988f) . which give time delays 
of a few x 10 s yr between the tidal interaction and the gas 
reaching the centre of the galaxy. 

This would be classed as a "minor merger", involving 
interaction of the galaxy with a dwarf satellite. Such interac- 
tions have been sho wn to be effective in feeding material into 
the nuclear regions (iHernauist fe Mihoil995h . This scenario 
for B1221— 423 would fit with the idea that low redshift radio 
galaxies are triggered by material falling into a pre-existing 
supermassive black hole, in contrast with very high redshift 
radio galaxies and quasars, which are often associated with 
major mergers and the simultaneous formation of a massive 
bulge and the central black hole. 

We thank Oliver Prouton for obtaining the optical im- 
ages of B1221— 423, and the anonymous referee for help- 
ful comments. RWH and EMS acknowledge support for this 
project from the Australian Research Council and the Uni- 
versity of Sydney Sesqui grant scheme. 
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